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Abstract 

We study an agent-based stock market model with heterogeneous agents and friction. 
Our model is based on that of [9] : the process of a stock price in a discrete-time framework 
is determined by temporary equilibria via agents' excess demand functions, and the 
diffusion approximation approach is applied to characterize the continuous-time limit (as 
transaction intervals shorten) as a solution of the corresponding stochastic differential 
equation (SDE). In this paper we further make the assumption that some of the agents 
are bound by either short sale constraints or budget constraints. Then we show that 
the continuous-time process of the stock price can be derived from a certain SDE with 
oblique reflection. Moreover we find that the short sale (respectively, budget) constraint 
causes overpricing (respectively, underpricing) . 

Keywords: Agent-based models. Liquidity problems. Short sale/budget constraints. 
Stochastic differential equations with oblique reflection. The Skorokhod problem. 

1 Introduction 

It is usual in mathematical finance to describe the price evolution of a risky asset such 
as a stock by a diffusion process. Geometric Brownian motion (GBM) is one of the most 
standard such processes for price fluctuation. Because of its simplicity and convenience, the 
GBM model is widely used in the context of option pricing/hedging, optimal investment, and 
many other financial problems. An important theme is to justify GBM from the economic 
viewpoint. For instance, a heuristic equilibrium argument for GBM is discussed in [21]. The 
justification of GBM as the rational expectations equilibrium is discussed in [2] and [16] . 

Recently there have been various studies of agent-based market models to explain the 
fluctuation of a price process. One representative study is the microeconomic approach of [8] 
and [5]: The process of the stock price is first given as a sequence of temporary price equilibria 
in a discrete-time market model with heterogeneous agents and then the price process in a 
continuous-time model is derived as the limit as the transaction time intervals shorten. Let us 
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introduce the outline of the model of [1]. Let / be the set of agents in the market and e^'*(p, u) 
be trader z's excess demand function for a proposed stock price p at the time k/n, where n G N 
and trades are executed at times t = 0, 1/n, 2/n, . . .. The parameter w is a sample point in 
the underlying probability space {Q,J^,P). The stock price process (5^/^)^=0,1,2,... is given as 
follows: At the initial time t = 0, the stock price is given by 5*0 = Sq. Then agents exhibit 
their excess demand e"Q(-,a;) and the next price is determined as a temporary equilibrium, 
that is, the solution p* of 

Y,eUp*{uj),u) = 0. 

After transactions at t = 0, the stock price changes to S']y„ = p*. Similarly, during the trading 
period k/n, the stock price before transactions is given by S'^/„ and agents' excess demands 

(e^'*(-, a;))jg/ make the stock price change to Finally, the process {S^j^k is given as 

the solution of 

E^^'(^(WH'^) = 0' fc = 0,1,2,..., S- = so. (1.1) 

For mathematical convenience we rewrite (11. ip by using the log price = log S*^ of the stock, 
obtaining 

E<'(^(WM'^) = 0' fc = l,2,..., Xo" = xo, (1.2) 

where e2'\x,uj) = e^'*(e^,Ci;) and Xq = logso. In [S] and the individual excess demand 
function e^'* is assumed to be given by 

i^^^) = (^.^) -^) +<^fc 

The parameter 5^'^{uj) denotes the liquidity demand, f^'^ix) is the reference level of agent i, 
and Fi is agent i's individual perception of the fundamental (log-)value. For a more precise 
economic interpretation of fll.3p . see Example 3) in Section 2 of [U]. In this case, fll.2p can be 
rewritten using the stochastic difference equation 

^r^+D/n - ^fc/n = mX;:/n - F) + 6^ /c = 0, 1, 2, . . . , (1.4) 

where 

En,i nn,i sr-~y n,i nn,i j-t sr-~y X^^^ 

i&I^k Pk ™ Z^iG/ "^fc Pk rn l^iel^k -n V^^">Vi i 

Pk= ' ^fc = , h = — , afc=2^afc (l + 7fe )• 

"fc "fc "fe 

Let (X,") t>o be an interpolated process of {X^^^)^ such as either 

Xr = X,"/„ te{k/n,{k + l)/n) (1.5) 

or 

Xr = {nt - A;)X(V,)/, + {k + 1 - nt)Xl,^ t E {k/n, [k + l)/n). (1.6) 
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Then, under some mathematical assumptions, the process converges to an Ornstein- 

Uhlenbeck (OU) process of the form 

dXt = /3{Xt - F)dt + adBt, Xq = xq. 

(This is a simphfied version of the result of [9]: They also treat a generalized process in a 
random environment.) This implies that the continuous-time stock price St = exp{Xt) is a 
geometric OU process. 

It is meaningful to consider the above diffusion approximation approach to derive the 
continuous-time process in a more general framework. A diffusion approximation for solutions 
of stochastic difference equation in the following form 

^(k+l)/n ~ ^k/n = ~]=^k{^k/ny ^) + ~Gk{^k/n-: ^)-: (l-^) 

\j n n 

with E[-F^(x)] = 0, is studied in [18] and [2l]-[26] under some mixing conditions. In [H], a 
case of the functional difference equation 

^rfc+l)/n - ^k/n = —f=Pk{{^r)r<k/n, ^) + "Cfc ( (X")^<fc/„, u) (1.8) 

is studied under strong mixing conditions and a certain additional dimensional condition. By 
using these results, we can apply the diffusion approximation approach to the agent-based 
market model for more general excess demand functions and derive several stock price models 
based on the framework of [9]. 

The aim of this paper is to construct an agent-based model of stock prices with market 
liquidity problems. In the real market, although there are agents who can buy and sell the 
stock freely to some extent, there also exist agents who cannot trade to their own satisfaction 
because of a shortage of cash, being prohibited from short selling, and so on. To consider how 
such a liquidity problem affects things, we construct a market model based on [9] under the 
following constraints: 

(I) Some of the agents cannot sell the more of the stock than the number of shares held (a 
short sale constraint), 

(II) Some of the agents cannot buy more of the stock higher than allowed by their budget, 
(a budget constraint). 

In each case (I)-(II), we will show that the continuous-time process of the stock price which is 
derived by shortening the transaction intervals is the solution of a certain stochastic differential 
equation with oblique reflection (SDER). Moreover, the value of the stock price under the short 
sale constraint is larger than it would be without such a constraint. The effect of the short sale 
constraint is discussed in [1], [5], [10] and the references therein, and our result is consistent 
with a common expectation, viz., that the short sale constraint causes overpricing (nonetheless, 
[1] pointed out that whether short sale constraints will always lead to overpricing is far from 
certain). On the other hand, we will also show that a budget constraint drives the stock price 
down. 

We now fix some notation. For an interval A C [0, oo), we denote by C{A] M"^) the set of 
continuous functions from A to and we use the abbreviations C| = C([0, T]; M'^) and = 
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C([0, oo); W^) (when d = 1, we simply write Ct and C.) For w G C^, we set \w\t = sup \w{t)\ 

0<t<T 

and |ty|oo = sup |if (t)|, where | ■ | is the Euchdean norm. We also define the following subspaces 

i>0 

of C^. 

C^,+ = {w = iw')UeC^; w\t)>0, 2 = 1,..., 4, 

^T,to = = i^li=i ^ ; ^XO) = 0' ^X^) is non-decreasing in t, z = 1, . . . , d}, 

and similarly and C^q for T = cxd. We introduce the canonical cx-algebra Bt = a{w{s) ; s < t) 
of C. We often consider the space of M^"''^-valued functions and then we start the index at 
zero, i.e., w{t) = {w\t))f^Q. The notation stands for the positive part of x, that is, 
= max{x, 0}. 

2 Main Results 

2.1 Model I: Stock Price Model with Short Sale Constraints 

Let / = {1, . . . , A^} be a finite set of agents who are active in the market which consists of 
a single stock. We assume that each agent i E I always stays in the market and no new agents 
enter. First we consider the discrete trading case, with market clearing times t = l/n,2/n, . . . 
for n G N = {1, 2, 3, . . .}. We denote by = {X^)t>o the log-price process of a stock, and by 
ip"''^ = {(Pt'^)t>o the number of shares of the stock held by agent i G /. The fiuctuations of the 
processes and cp"' = {ip'^'^)f^^ are found as follows. At the initial time t = 0, every agent 
i G / has (Pq'^ = > shares of a stock and the initial log-price of the stock is Xq = xq G M. 
After trading is finished for the period t = k/n, the log-price X^, t G {k/n,{k + l)/n] is 
determined by the market clearing condition 

Y.'riXi,^l)/n,X-,U-^lfj=0 (2.1) 

and linear interpolation (11. 6p . where e^'^{x, w,u]ip) : MxCxfi" xM — > M is the excess demand 
function of i G / at t = /c/ra on the underlying probability space (f2", J-'"', P"'), and ip^'^ is the 
quantity of stock held by i. Here e2'\x,w,u] ip) is assumed to be B^R) B^/n ® J^"' ® Im- 
measurable, so that 

el'\x,w,u;(p) = e5!'*(x,M;(- A {k/n)),u;(p). 

After determining the log-price up to t = (/c + l)/n by (12. ip . the agents' holdings of the stock 
are set by 

f{k+l)/n = fk/n + Sfc (^(fc+l)/„, ^ , ^, Vk/n) (2-2) 

and 

= {nt-k)ip'l^^-^.,+{k + l-nt)ipl']^, t e {k/n, {k + l)/n) (2.3) 
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ioT i = 1, . . . , N . We remark that fl2.ip implies that the total number of shares of the stock in 

N N 

the market is constant, i.e., = for any n and t. 

i=l i=l 

All agents have their individual excess demand e2'^{x,w,u) before considering friction, 
but they do not always exhibit e^'^{x,w,u) itself to the market. Agents are divided into two 
groups, I = hU h, where Ji = {1, . . . , iVij and /2 = {A^i + 1, . . . , iV} for some 1 < A^i < iV. 
Agents in the first group are prohibited from selling short, thus their excess demand aXt = k/n 
is not lower than v^^^*„ and the process of their stock holdings is always non-negative. Agents 

in the second group I2 are allowed to sell short, so they can exhibit the value e^'*(-, X", w) 
itself as their excess demand. Then the excess demand function is finally defined as 

^n,i/__ f max{e^'*(x,ty,w), -v?}, i e h 



el\x,w,uj-^) = {Z7( '\ '^' (2.4) 



We divide e^'*(x,w) into two parts. 



Here 



e^'^fx, w) = fk'\x, w) + g^'^iw). (2.5) 



fk\x, w) = el'\x, w) - el'\w{k/n),w), gl'\w) = el\w{k/n),w). 

The function g^'^{w) is the excess demand when the price of the stock does not change. We 
further divide g]^'\w) into 

9k M = -gk M + -^gk (^)> 



where g^^ivS) is a mean zero random function. We remark that g^^ ^ g^'^ and g^'^ are all 
-Bfc/n ® -^"-measurable. The function f^'\x, w) is an additional excess demand associated with 
the change of the price. Usually the higher the stock price becomes, the lower is an agent's 
demand for the stock. Thus it is natural that fj!'\x,w) is supposed to be decreasing in x. 
Furthermore, we assume 

[Al] The function f^'^{x,w) is deterministic (i.e., independent of a; G fi"), continuous in 
{x,w), and three times continuously differentiable in x. Moreover, there exist positive 
constants Kq and 60 such that 

-Ko<^f:^'\x,w)<-6o. (2.6) 

Furthermore, g^'^{w) and g^'^{w) are continuous in w almost surely. 
Condition [Al] implies that when I2 is not empty (that is, when A'^i < A^), the equation 

J]er(a:,i.;^^) = (2.7) 
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has a unique solution x for any fixed w and (v?*)^, since the left-hand side of fl2.7p is strictly 
decreasing in x. In fact, the existence of a unique solution of (12 .Tp is also guaranteed even if 



N 



I = Ii, provided > 0. However, hereafter we always assume Ni < N for some technical 



1=1 



reasons which we will explain later. 

By (12.11) and (12.21) . we can construct the processes {X^)t and {(Pt)t- We are interested in 
the limit of the (1 + A^)-dimensional process = (X", (pt'^, . . . , ^Pt'^)t as n — )■ oo. More 

precisely, we consider the weak limit of the distribution /i" = P(H" G ■) on C^~^^ . 

We will define more conditions. 

[A2] For every M > 0, there exists a positive constant Cm > such that 



3 Ql 



^ sup 

;_0 \x\,\-w\oo<M 



+E"[ sup i^rHi'i+E"[ sup i^rHr']<CM 



\w\^<M \w\ac<M 



for any n,k & N and i = 1, . . . , N, where E" is the expectation with respect to P" (we 
simply denote this by E when there is no possibility of confusion). 

[A3] The cr-algebras cr{g^'\ g^'^ ; i = 1, . . . , N), k = 1,2, . . ., are independent. 

[A4] Let 

ifH = ^,r,\w{kln),w), aY\w) = E"[^r(«^)^r^(^)]- 

For each i,j & I the following limits exist 

f3\t,w) = lim l3^;Uw), i\t,w) = lim T^H, 
a'^{t,w) = lim ay^Jw) 

r— >-oo ^ ' 

uniformly on any compact subset of C and for any t > 0, and 

a'(t,w) = hm «|^J](m;) 

uniformly on any compact subset of [0, oo) x C. 
[A5] Define Q^iw) = {QY^ {w))i<i<M,i<j<Ni as 

QV\'Uj) = (1 - 5ij)al'\w)al\w), 
where 5ij is the Kronecker delta and 

TV 

«^(^) = E<'H, = l/{al{w)-al'\w)). 

i=l 

There exists V = {V'^)^j=i E R^' ® M^^ such that V'' = 0, Ql''\w) < V'^ for each 
i, j = 1, . . . , Ni and p{V) < 1, where p{V) denotes the spectral radius of V. 
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Note that the inequahty A''i < is essential for condition [A5]. Indeed, if A''i = A^, then 
the calculation 

N N N N 

1=1 lm-l=i ll = l * = 1 

N N 

'm-l=l h = l 

indicates p{QkH) = lim lHQkHrWl^"" = 1, where {{QkHrY^ is the z,jth element of 

m—^oo 

N 

the mth power of the matrix Q2{w) and the norm || ■ ||i stands for ||A||i = maxN for 

1=1 

A = {A'^)ij e (g) M^, thus we cannot find such a matrix V in [A5]. 

[A6] Let a{t,w) = {a^\t,w))^j^-^ be an A^-dimensional matrix-valued function such that 

N 

a^^{t,w) = a^'^{t, w)a^"^{t, w). For any T > 0, there exists a positive constant Ct 

m=l 

such that 

\^\t,w)\ + \fit,w)\ + \a'\t,w)\ < Cril + \w\t) 
for each i,j = l,...,N,0<t<T and w G M, where 

p\t, w) = f3^{t, w) + fit, w), fit, w) = 7 ; ,2 E 

za[t,w) ^^^^^ 

We now introduce an SDER. 

N N-L 

dXt = lPit,X)dt + ^a^^it,X)dBi + ^a^it,X)dL{, Xo = xo, (2.8) 

N Ni 

d^\ = h\t, X)dt + E ^'Ki^ X)dBi + li,ii)dLi - E Q'^it, X)dL{, ^i = ^\ i = l,...,N, 

j=i i=i 

where 

Q!*(t, w) = w) - Q!*(t, w)), Q^^it,w) = (1 - 6ij)a\t,w)a^it,w) 
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and U{t,w), a^^{t,w), < i < N,l < j < N are given by 

Y,^'{t,w)/a{t,w) (2 = 0) 

TV 

P\t,w) -a'{t,w)^/3^{t,w)/a{t,w) {i > 1), 



(2.9) 



i=i 



N 



fc=i 

TV 

a'^{t, w) - a\t, w) J2 ^^^ (^' w)/a{t, w) {i > 1) 



(2.10) 



k=l 



N 



with a(t,w) = ^^a*(t,w). We say that an (1 + + A?"!)- dimensional continuous adapted 



2=1 



stochastic process (Xf, (pt, Lt)t = {Xt, {^pD^i, (LD^^i) is a solution of fl2.8p on a given filtered 
space (fi, J-", {J^t)t,P) equipped with an A^- dimensional (J-'Ji-Brownian motion Bf = {BD^^^ if 

. P(^GC^)=P(LGCfo^) = l, 



^ {ipi.>o}dLl. = ioT i = 1, . . . , Ni almost surely, 
The processes {Xt)t, i(ft)t, {Lt)t and {Bt)t satisfy 



Xo+ f h\r,X)dr + y^ f a'''J{r,X)dBi + y^ f a'{r,X)dLi, (2.11 
Jo Jo Jo 

/ b\r,X)dr-y2 / a'\r,X)dBl + h,{{}Li-y2 / Q''{r,X)dLi 
Jo -^^ Jo Jo 



for t > and z = 1 , . . . , almost surely. 

We call the process L = {L\)i(zi^^t>o a regulator associated with S = {X,Lp). Now we present 
our final assumption. 

[A7] A solution of f l2.8p is unique in law. 

For instance, if dfj!/dx, i = 1, . . . ,N, is constant (and so is the matrix Q), condition [A7] 
holds under a Lipschitz condition on the coefficients 6* and a*-' (see [1] for instance: although 
the form of our SDE (12.81) is a little special, the arguments in these papers also works.) For 
other sufficient conditions for [A7], see [7] and 

By [A7], the distributions /i = P(E G •) and u = -P((S, L) G ■) are unique if exist. We are 
now prepared to state our main result. 

Theorem 1. Assume [A1]-[A7]. Then /x and u exist and the distribution /x" converges weakly 
to /i on C^~^^ as n ^ oo. 
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The proof is in Section [31 Theorem [T]imphes that the hmit of H" = (X", v^"'^, . . . , f/)"' ) 
is characterized as the solution of an SDER, viz., (12. 8p . The regulator process L prevents the 
shares of stock of an agent i from taking a negative value. The infinitesimal term li^{i)dL\ in 
(12.81) works only when the agent i E h hopes to sell more of the stock than they hold. 

Here we consider the case where all agents in the market can sell short (i.e., Ni = 0). This 
is a version of Theorem 1 of [H] and converges weakly to the unique solution X of 

JV 

dXt = h\t, X)dt + ^ (T°^ (t, X)dB'l Xo = xq. (2.12) 
j=i 

In this case, no agents are bound by the short sale prohibition, and the process {Xt)t represents 
the log-price of the stock without friction. 

On the other hand, when Ii is not empty, agents in Ii may not be able to exhibit their 
primary excess demand (which is described as e^'* in the discrete-time model.) The (log- 
)price Xt is pushed up by the gap between the actual excess demand with the primary excess 
demand, so Xt is larger than Xf. To describe such a phenomenon, we give the following 
additional assumption. 

[A8] There exist measurable functions a^{t,x), j = 1, . . . , A^, such that o"°"'(t, w) = {t, w{t)) 
for each t and 

JV 

Y,\^'it,x)-a'it,y)\<pi\x-y\), t>0, x,y gM^ 
i=i 

for some strictly increasing continuous function p : [0, oo) — )■ [0, oo) satisfying p(0) = 
and 

>/(0,oo) P(0^ 

[A9] The pathwise uniqueness of solutions of (I2.12p holds. 

Theorem 2. Assume [A1]-[A9]. Let {X,ip,L) (resp., X) be a solution of h2. S\) (resp., h2.12^) ) 
on a given filtered space (f2, J-", {J^t)ti P) equipped with an N -dimensional Brownian motion B. 
Then Xt > Xt for any t > almost surely. 

Theorem [2] is directly obtained by the same proof as Theorem 1.1 in [12]. This suggests 
the assertion that the short sale constraint causes the overpricing in our model. 



2.2 Model II: Stock Price Model under Budget Constraint 

We also consider the case where some of the agents, Ji = {1, . . . , A'^i}, cannot borrow cash. 
In the previous section, an excess demand function with no friction e^'*(x,ty) is understood 
as shares of the stock which an agent wants to buy. In this section, we interpret e^'*(a;,w) to 
mean an excess demand in terms of dollars. We also define 



eY{x,w]y) 



mm.{eY{x,w),y}, i e h 
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The market clearing condition is now expressed by 

N 



(2.13) 



where is the agent's amount of cash held at time t. Then the process of the log-price of 
the stock {X^)t and the cash holdings of an agent i are given by f ll.6p . fl2.13p . 



W, 



{k+l)/n k/n 



k/n> 



(2.14) 



and 



W-^' = {nt - k)W;il,y^ + {k + l- nt)W-ll t G {k/n, {k + l)/n). (2.15) 



We assume that e^'* has the same form as (12. 5p . We also assume [A1]-[A7], replacing b, a, 
and (123) with 



b\t,w) 



Y,PKtM/a{t,w) (z = o) 



N 



I3\t, w) + a\t, w) E /3^'(t, w)/a{t, w) (i > 1), 

^a''^{t,w)/a{t,w) (i = 0) 

=1 

TV 

a'^{t,w) + a\t,w)^a''^{t,w)/a{t,w) {i > 1) 



N 



k=l 



k=l 



and 



N Ni 

dXt = b\t,X)dt + J2^^^it^X)dBi-Y^a^{t,X)dLi, Xo = xq, 

TV TVi 

dWi = U{t, X)dt + E a'^it, X)dBi + 1/, (z)rfLj + ^ Q'^it, X)dL{, = d 

i=i j=i 



(2.16) 



z = l,...,iV, 



where c* > is the initial cash holdings of agent i. Then we have the following theorem. 

Theorem 3. The distribution of'E"' = (X", W"''^, . . . , H^"'^) converges weakly to a solution 
of( ^rm> on C^+^ asn^oo. 

Theorem 4. Let E = {X, W^, W^) (resp., X) be a solution of (resp., ^2JE ) on 

a given filtered space {Q, J^, {J^t)t, P) equipped with an N- dimensional Brownian motion B. 
Then Xt < Xt for any t > almost surely. 



We omit the proofs of Theorems [3H1] since they are almost the same as those of Theorems 



mm 
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3 Proof of Theorem [T] 

Take any M > |xo| and let e C°°(M; [0, 1]) be such that z/^a/ (y) = 1 on |y| < M/2 and 
i^Miy) = on \y\ > M. We set 

e1'^''\x, w) = -(1 - ■ipM{'w{k/n)))a1'\w){x - w{k/n)) + ipM{'w{k/n))el''{x, w). 

We define X'^'*'^ (^"'^^'* and el^^'\x,w;^) similarly to O-Q, replacing g^'^ with el^'\ 
First we consider the convergence of the truncated processes S"'*^ = (X"'^, (/j"'^^'^, . . . , y?"'^'^), 
n G N for fixed M. We can easily see the following proposition (Proposition 1 in [T^). 

Proposition 1. For any lo, if < M, then |X"'*^(a;)| < M for all r e [0,t]. 

We rearrange our market clearing equation into the form of a difference equation. Since it 
follows that 

c~n,M,i/ \ -| ~n,M.it \ , i ~n,M,i/ \ \ 

maxje^t' ' (x, w), = e^' {x,w) + [-e,: + , 

we get 

N N Ni 

i=l i=l i=l 

where = (— e^'^^'*(X('^';Jf^^y^, X"'^^) — (^9^^^^'*)+. Using Taylor's theorem, we have 



^(fc+l)/n ^k/n 



1 ^ / 1 



EC'''' (3.1) 



i=l 



and 



n,M,i _ n,M,i 
^{k+l)/n ~ ^k/n 



/ vn',M\ f n,if-crn,M\ , ^ ^ n,i / ■c^n,M\ / vn,M vn,M\2 , n,M,i 

- + hA^ftl''^ ^ (3.2) 



where 



n,M,i -'- /" /I \2 ^ £n,i/ x^n.Af , /i \ ^n,M \rn.M\ j / \rn,M ■srn,M\3 

=2 (1 - d^fk («^{fc+l)/n + (1 - ^)^fc/n ' )^«(^{fc+l)/n " ^fc/n ) ' 

Substituting f l3.ip into itself and into fl3.2p . we get 

1 f ^ iVi >| 

^(fe+l)/n ^fe/n - ^nf-^n,M\\Z^^k \^ I^Z^'h p 



n,A/,i ,n,M,i _ TTn,M,i ( vnM\ "fc l"^ J ] un,M,j I vn,M\ i \^ ^'^.^J 



+^iMVk 



Ik 1 
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where 



1 1 



n.i/ \ N 

' j,m=l 

~n,M,i _ I / ■^n,M\ I n,M,i . Ik \^ J -n,M 



I f yn,M\ I n,M,i I Ik \ I ^n, 

Vk jSl^Vk +2V^m(X,;J^ ( -=(7^' (X ' )+p,' 

i=l ) L i=l i=l ^"^ 

I V jj^;^ \i = l 



n,M,i 
Pk 



n,M,i ^iT't^ f vn,M\ , -'- , n,i / -v^n,!;? \ / -vAn,M -vAn,M\2 , n,M,i 

Pk - -9k ) + 2^k )y^ik+l)/n- ^k/n > +^k 



n 



Thus, if we set ?7^'^'* = (1 — a2'\w) / a'^{w))ri^ 



k 5 



[nt]-l [nt]-l 
fc=0 fc=0 



and 



N 



= -0 + E /„ -^wk^^^r'^'N (3.3) 



T^n,M,i yn,M,i \^ / L^V^MJ ^- 1 Ar /'Q/l^ 

^0 '^[nr]l^ ) 



then 



-n,M,i 



(3.5) 



1=1 

The equahty (jM]) seems to imply that ((/?"'^-^'\ L"'*^'*)^\ is a solution of the Skorokhod 
problem with oblique reflection in the non-negative orthant associated with (y"'*^'*)^\ (see [1], 

POO 

[11], [19], [20] and [22]). However this is not strictly true, since the equality / ip'^'^^^'dL'^'' = 

Jo 

does not hold by virtue of the linear interpolation (12. 3p . The following proposition tells us that 
('/'fc/^'*' -^fc/^'*)£i' ^ ^ were, a solution of the corresponding Skorokhod problem 



k/n ' fc/' 

in discrete-time 
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Proposition 2. For every k G Z+, 



^^^'' = o^f<i.UZ/ • (3-7) 

Proof. It is obvious tliat tlie left-liand side of (13. 7p is not less tlian the right-hand side. We 
suppose 

/ ^1 .l/n \ 
^Ifn' > [T. I ^?F4(^"'')^^"*'' - "^Sn'") (3-8) 

for / = 0. . . . , /c. By (13. 8 p with I = k, we have 

pk/n 


This inequality gives = (by the definition of ffl'^['^) and thus L^j^'^iy^ ~ ^ifn^- Using 

(13.81) again with / = A; — 1, we similarly get rj^'^'^ = 0. Inductively we see that = Q for 

t G [0,k/n] and this contradicts (13. 8p . Then we obtain the assertion. ■ 

Using the above proposition and the same arguments as Theorem 2 in [22], we get the 
following proposition. 

Proposition 3. For every < I < k, 

Ni Ni 



y \Llf'' - L"/"'f < max y |F"f - K^'^'f 

/ J I k/n l/n I — l<rn<k ' ' "^Z" " 

i=l ~ ~ i=\ 



for some K > depending only on V . 
The equality (13. 2p also indicates 

n,M,i ■(yn,M,i , fn,MA ■ i at /o r\\ 

^k/n =^k + Lk/n ' , ^ = 1, . . . , A^i, (3.9) 

where 

= 2^\VJM[^i/n ) [9l >^2^l )y^{l+l)/n- ^l/n > + 

1=0 { ^ 

-ar(x-^)(xjvV-^r/f)|. 

fc=0 

The equality (13. 9p corresponds to the classical Skorokhod problem for each i = l,...,Ni. 
Similarly to Proposition |2l we obtain the following. 
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Proposition 4. For every k G L'^'j^'^ = (~^//n^^'*)+- 

Next we evaluate the moment of the process. 
Proposition 5. EflXj^'f^^/^ - Xlf^\'^^] < Cm/u^^ k e Z+ for some Cm > 0. 
Proof. Set 

N 

7-f4.f\^n.M •trn.M\ I \^n,AI\ 7/ \ \ ^ n,M,i/ \^n,M n.M,i\ 

fit) = /(^(^(fc+l)/n - ^k/n ) + ^k/n )' f i^) = l^^k ' (a^, ^ ' ; V^fc/^ )• 

i=l 

Then we have /(I) = and 

A^i N 



Using the fundamental theorem of calculus, we get 

/(o) = /(o) - /(I) = - r /'(t(x;'^i)/„ - x-f ) + x-f - x-f ). (3.10) 

t/ 

Since > for i = 1, . . . , Xi and /'(x) < —N25q, we get 

i^5f,/„ - nfj <W<^X (3-11) 

By [A2] and Proposition [H we have 



E[^M(X,";f )|^,"'^(X"'*Ori = E[ sup \gl'\y)r] < %. (3.12) 



Our assertion follows immediately by (13. lip and ( I3.12p . ■ 
The above proposition and [Al]-[A2] lead us to the following. 

roposition 6. Y\\e^: ' (X^^'^^^/^, X • )| ] + E[|%' 'I ] + E[|%' 'I ] + Efk^' ' T] < 
Cu/n^'^ , k G Z+, /or some Ca/ > 0. 

Proposition 7. EflL^^'^^l*^] + Eil^)"'*^!^] < Ca/,* /or some CM,t > 0. 

Proof. It suffices to estimate E[| -^j^t]'/^, H • By [Al], Proposition HJ and the equality 

i=l 



^{k+l)/n ~ ^k/n - ~ ^k/n > ^ {\X;^f^ \<M} I'^-'^'^J 
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(obtained by the definitions of el'^^'\ fi^'^'\ and X"'^''^), we have 



EL?'^/^ < max 
[nt]/n - 



i=l 



i=l 



rin,M,i r^n,M,i . ^ ST^ I 7n,M,i ) 7 



1=0 ' ^ 



< 



max \GJ + max G, ' 
Vo<fc<H " o<fe<H " 



N [nt]-l 



M,j| 



i=l fc=0 



fc=0 



where 



n,M,i 



G 



n,M,i 



V , = 1 ;=o ' ^ ^ 



' n,M,i 



(fc+l)/n ~ "^fc/n ^ + ^^k 



Since 



< 



by virtue of [Al], we get 

Ni Ni 



i=l ^ 



i=l 



max G,, 

0</c<[nt] 



fc 



i=l fc=0 



Efn,M,i 



i=l 



and thus 

Ni 

y ^ff 

/ J \nt\/n 



i=l 



< 



V ^2^0 J I \0<k<lnt] 



max I G 

0<fc<[nt] 



n 



N [nt]-l 

E E K 

i=l k=0 



By [A4] and Propositions ISHHl we have 



N 



i=l 



[nt]-l 

Ei'^: 



n,M,i I 



k=0 



< G 
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for some C'j^j^ > 0. Moreover, since and (G^'*^'*)^ are both (^^)fc-martingales, where 

QJ^ is the (J- algebra generated by (7"'* and for z = 1, . . . , A^, / = 0, . . . ,k — 1 (note that Qq 
is a trivial a-algebra), the Doob inequality implies 



E[ max |G^'*-''f + max IG^'^'^'T] < 

" " " 0<k<[nt] 
nt]-l 



n,M,i 1 81 



< 



0<k<[nt] 

2 /8 



n2 V7 



8-1 



< r" 



fc=0 



(3.16) 



for some C'lj^ > 0. Now by fl3.14l) - fl3.16p . we obtain the assertion. 
Proposition 8. For every t > there exists a constant CM,t > such that 

[nt]-l 



~n,A/,j 



< 



c 



M,t 



k=0 



n 



Proof. It suffices to show this for e^'*^ instead of e^'^^'*. A straightforward calculation gives 



[nt]^l Ni [nt]-l 



[n,f]-l 



4i 



fc=0 



< 



max 

0<fc<[nt]- 



fc=0 



E Ei 



n,M,i 
[nt] 



1=1 



max \wj 
o<fc<H-i 



[nt]-l 



M|4i 



fc=0 



where 



7n,M 



n.M 



Ni ^ / 1 

5:c"'^+2^M(x-f)$:(^r(^ 



J Pk 



N 



N 



n,M\2 
k/n ) 



^ j,7 = l \i=l 



M,: 



rn,M\ 81 



< 



Since [A2] and Proposition imply E[V^a/(^"' )(pfc' '')^] < C'j,j/n^ and thus E[(^; 
^m/''^^; E[(^fc'^)''] — ^m/''^^ some Cj^^C'l^ > (by virtue of Proposition E]), using 
Proposition [TJ we get 



[nt]-l 



fc=0 



it] 



i=l 

1 1 
= + - 

n n 



1/4 



[nt]-l 

E E[l*' 



1/4 



A/ 1 81 



A:=0 



for some Cj^,f^,Cj^jf > 0. This implies the assertion. 
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Proposition 9. E[ sup \Z^'^'\^] < Cm,t for some Cm.t > 0. 

0<t<T 

Proof. Our assertion is obtained from [A2], Proposition [HI fl3.16p . by calculating 



E[ sup 

0<t<T 



< E 



max 

0<fc<[n.T] 



n 



[nT] 



< C'MT{n max \Gl'''f] 



~n,M,i 1 4l 



n 



;j=0 l«'l=o<M 



A:=0 



for some C'j^.jrp > 0. ■ 
Proposition 10. Let (^fc)fc be uniformly bounded random variables such that is G^-i' 



adapted. (Here, {Gk)k is defined as in the proof of Proposition^) Put Z. 
Then {Z^'^^)t is tight on C. 
Proof. It suffices to show 



■n,M,i 



j.pn,Af,. _ ^n,M,i|2|^n,M,i _ < ^^^^^^^ _ ^^3/2^ 



and 



limsupEfl^"'""'" - ^r""'? ] < CmA^ - 



(3.17) 



(3.18) 



for any < s < t < n < T and for some Cm,t > (cf. [I5]). Set $ = - Z^'^^''\ for 

brevity. Using the inequality 



(X^^O = X^^? + 2X^a;i(xi H hx;), xi,...,XfcG 

Z=l 1=1 1=1 



and the uniform boundedness of (^fc)fc, we get 



< 3E 



[nu] — 1 



^k=[nt] 
[nu] — 1 



(3.19) 



< C'{ (E[(i/p4'r]^/^ + E[(i/jf^)Y/^+ E E[(i/r'?]^/M E[$T/^ 



fe=[n.i] 



[mi] — 1 



I \ ^ l-r^l^n Tjn,M,i/ -trn,M\ f ryn,M,i ^n,M,i\^-\ 

+ |E[4fc^fe (A )(^fc/„ -^[„t]/„)^J 

k=[nt] 



(3.20) 
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for some positive constant C . By [A2] and Proposition |H1 we see that 



n 



(3.21) 



for some C'lij, > 0. Moreover, [A3] implies 



and hence 



A 



< E 



-\pM[,J^k/n J'^k )+^k 



1/4 



E[ sup v-fe/. 

0<fe<[n.u]-l 



< 



-^r4/0T/'E[$T/^ 

CZt (- + E[{e2'''nV) E[ sup |Z,"'^'^'f e[<I>2]V', ^ < M - 1 (3.22) 

' \n J 0<r<T 



for some C'J^j. > 0. The inequahties fl3.20p -( 13.22p and Propositions [SHni imply 



[nu] — [nt] + 1 



n 



E[$ 



2ll/2 



for some C^/V > 0. Replacing $ with 1 and performing the same calculation, we have 



E[|Z 



n,M,i yn,M,i\2-\ ^ r^llll 



[nt] - [ns] + 1 
n 



hence 



[nu] — [ns] + 1 \ '^^ 



n 



(3.23) 



(3.24) 



The inequality f l3.23p immediately leads to fl3.18p . The inequality (I3.17P is now obtained by 
fl3.24p and the same argument as in the proof of Theorem 14.1 in [3]. ■ 

As a consequence of Propositions [3] and [TOl we can see the tightness of the processes with 
fixed M. 

Proposition 11. A family of processes (X"'^, y"'^^ Z"'^^ L"'*^ L"'*'^)„ is tight on C^Ci+^+^O, 

Proof. The tightness of (y"-'*^^ i^n^M-^ obtained directly from Proposition [TOl Then Theorem 
7.3 in [3] implies 



lim lim sup P{wt{5; y"'^'*) > e') = 0, i = 0, . . . , iV 

(S— >0 n— ^oo 



(3.25) 



for every > and T > 0, where wt{5\ x) stands for a modulus of continuity, i.e., wt{5\ x) 
sup |x(t) — x{s)\. 

Q<s<t<T,\t-s\<S 
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For a while, we set n sufficiently large so that 1/n < 6. Let 0<s<t<Tbe such that 
\t- s\< 6. By [Al], we get 

lfn,M,i _ fn,M,i\ ^ -^n,MA , f,n,M,i , \^ (ij "fc ) 1 

■•n,M,i , / 1 , -^0 \ „.. /r. Tn,M,i\ 



< 2 max C 1+ ,,, \wT{6;L-'''n- (3-26) 

0<fc<[nT] V (iV - l)5o / 

Similarly, Proposition |3] implies 

iVi 

|^n,M,i _ j^n,M,i\ < 3 max 7/"'^^'^ + vIYwt(S; F"'^^'^). (3.27) 
By fl3.27p . Proposition El and the Chebyshev inequality, it follows that 

max > £^/2) 



j=l ^ V U 1/ ^^Q 

< yP wT{S- A 



M,j\4i 



256CM([rir] + l) 



for any e > 0. Taking limsup, letting 5 — 0, and applying f l3.25p . we get 



lim lim sup P{wt{6] L"'*'-*) > e) = 0. (3.28) 



Similarly, fl3.27p . ( I3.28p . and Proposition [6] imply 



lim lim sup P{wt{5] L""'^^'') > e) = (3.29) 



t s \ — \ t s I (M_^\X^Z-^' * « I 



for any e > 0. Furthermore, the inequality 

Ni 

gives 

lim lim sup P{wt{S; X"'^) > e) = 0. (3.30) 

Our assertion is obtained from f l3.28p -f l3.30p and the fact that the initial values Xq'*^, Lq'^^'*, 
^n,M,i constants. ■ 
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PropositiondHtells us that for any sequence (nfc)^ C N there is a subsequence (rifcj/ C {nk)k 
such that F^^r*^, Z"'';'*^ L"'';'^'^ L"'=i'^'^) converges weakly to a certain continuous pro- 

cess (X*^, Z*^ L*^ L^O defined on some probabihty space (fi^^J^^^P^O on C'^+^+^O 
as / — oo. Furthermore, since [Al] and [A4] imply that the convergences 

are uniform in t G [0,r] for all T > 0, using Proposition [71 Theorem 2.2 in [17], and the 
continuous mapping theorem (Theorem 2.7 in [3]), we obtain the weak convergence of 

^^UkjM^ ^n'k,,M^ Y^ki'M^ Z^k,,M^ JJik^M ^ frikiM ^ jnkiM'^ 

(X*^ (^*^ Y"', Z^', J^) , / ^ oo, (3.31) 

where 

fn,M,i f ~n,i / •srn,M\ j t n,M A Jn,M,i f /^n.ij / ^n,M\ j t n,MA 


Jo Jo 

and 

Ni 

V't +1/1 W^t Z^^t • (3-32) 

Note that ([S3]) and f[X^ tell us 

= + Jf '\ (3.33) 

Let us introduce a filtration on (fi*^J^^^P^O• We define = (T(Xj^^ , Lf ; r < t) and 
let ( Jf^)t be an enlarged filtration of (^f^)* such that T^' , (J^t^h P^^) satisfies the usual 
condition. We notice that the processes (/? , and L^'^ are also (J-'/'^)radapted. Now let us 
define 



^0 

where p'''\t, w) = iljM{w{t))P\t, w) + tljM{w{t))^\t, w). 

Proposition 12. For each i, j = 1, . . . , N , the processes {Nf^'^)t and (iV/^'*-')^ are both (J-'/^)f- 
martingales. 

Proof. It suffices to show that [N^ ^'^)t and (iV*^'*'')^ are (^j*'^)rniartingales. Set 

fc=0 

/z^'X^) = E[/.^'^'^(^)], /^r'^'X^) = i^M{w{k/n))hr''\w). 
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Then we have 



E 



sup 

0<t<T 



where 



k-l 



rn,A/,j 



max ( H, 

0<fc<[nr] + l 



[nt] 
k=0 



Tjn,M,i \ ^ , f vn,M\/ ,n,M,i/ ^n,M\ r,n,M,i/ yn,M\\ 



(=0 



Then the same calculation as fl3.16p and Proposition [8] leads us to 



E 



sup 

0<t<T 



0, n — 7- oo. 



(3.35) 



Moreover [A4] implies lim '^{w) = (3 '\t,w) uniformly on any compact subset of C for 
all t > 0. Thus, using f l3.3ip . we get 



I — )■ oo 



(3.36) 



weakly on C. The convergences fl3.35p - fl3.36p and the continuous mapping theorem imply that 
^^n,,,M^^n,,,M^^„,,,,M^^n,,,M) ^ (X^^ Z*^ L*^ iV^^) , / ^ OO Weakly ou C^^. (3.37) 

On the other hand, by [A3]-[A4], we have 



and 



^* E[^M(X,"'*0'«H(^"'*0*]^r 



{nt - [nt]){[nt] + 1 - nt) 
n 

{ns — [ns]){[ns] + 1 — ns) 



E[V^M(Xp4;ja-f(X-^) 



$1 



E[V^M(X("4;)aH(X-^0$] 



n 
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for any mGN, 0<Si<---<s<t and any bounded continous function $ : — y 
M. Thus, using [A2], fl3.37p . and the dominated convergence theorem, we obtain 

E"^[(iVf - N^'nmX^X=i^ (Zif iL^X=i)] = 0' (3-38) 
= E^^[ /*^M(Xf )^a'^-(r,X*Orfr$((Xi;0;^i, (^?)r=i, (3-39) 

J s 

(where E*''^ stands for the expectation under P*''^), which imply our assertion. ■ 

By Proposition [T2l and the martingale representation theorem (Theorem 3.4.2 in [I3]), we 
can construct an enlarged filtered space {Cl^",T^^, {^l^)t,P^^) of (fi^^ J^*^ (j;^)f,P^) and 
find an iV-dimensional (J^/^^)t-Brownian motion {Bf^)t such that 

= E / ^M{X^'V^{r,X^)dB^'\ (3.40) 

where the stochastic processes on (fi^^, J^^\ {J-'l^)t, P^'^) are regarded as defined on (fi*^, J^^^, 
{P^^)t, P^^) canonically. Moreover the process {Z^^'^)t becomes an (^/'^)j-semimartingale and 
we can define the stochastic integral 

Jo Jo ^-^^ ^0 

for an (jf ) (-progressively measurable process {C,t)t (under suitable moment conditions). 
Proposition 13. The following equalities hold. 



i=i 



Proof. By Proposition [S] and 

N [nt]-l 

nTI ; / vri.M\ ~n,i / ^n.M\\2 , I / / vn,M \ , n,M.i / x^n.M\ \l ^ ^ 



iV [Utj — 1 

-EE E[i^M(x-f)^-^(x-^op+ < 



n 

i=l k=0 



for some Cm,* > 0, we can apply Theorem 2.2 in [T7] to arrive at the weak convergence of 
Our assertion is now obtained using fl3.3p - fl3.4l) . fl3.4ip . and the continuous mapping theorem. 
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By fl3.31l) - fl3.34l) . fl3.40p . and Proposition [T^HT^ we obtain the following proposition. 
Proposition 14. The pair [X^'^ ^tp^) (and the regulator process L^^ ) is the solution of 

N Ni 

dX^ = b'''\t,X^')dt + Y,^^''''it,X'')dBl''' + Y,<!^'it,X^')dLf^', X^ = xo, 

N 

d^f'' = b'''\t,X^)dt + Y,&^'''it,X''')dB^''^ + li,{t)dLf'' 

-J2Q''{t,X^)dLf''\ (^^ = $\ z = l,...,N, (3.42) 

i=i 

where o ^ and a^^ are given by (12. 9p - (I2.10p upon replacing /3* and a^^ by l3^^'^ and cr*^'*-' . 

Proof. It is obvious that {ipf^'^)t is non-negative, (Lf^'-')t is non-decreasing, and L^'-' = for 
i = 1, . . . , N and j = 1, . . . ,Ni. The rest of the proof is to show 

oo 

i^f^'dLf'' = 0, « = 1, . . . , A^i almost surely. (3.43) 







By the definition of L"' , we have 

J, [nT]-l 

n,M,i irn.MA \ ^ n,M,i n,M,i , / t-i r rrl\ n,M,i n,M,i T^\n/0/i/i\ 

^i[nr]+l)/ndW = ^ ^iUl)/nVl + (^^ " [^^])V5([nr]+l)/n^[nT] =0' T > 0.(3.44) 

Propositions EHI] imply 

t 

n,M,i ,„n,M,i\^Tn,M,i\ 



E[ sup I [\^l%y^ - ^-^'0^^: 
0<t<T Jo 



N Ni 

< EP^'*''*)T/MFe[ max (i/,"'^'^y]V' + 'rE[ max (C''''^')^/' 

AT /[nT] \ hnT\ ^ '''' 

j=l \k=0 J j=l \k=0 
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C A f T 

< — ^ — ^ 0, n ^ cx), T > (3.45) 

for some > 0. Using ([331D-(l335D and Theorem 2.2 in [IT], we obtain <^M)- ■ 

Here, (E33) and (EH imply that ((</?^^'0£\, (iv^'0,=\) is a solution of the Skorokhod 
problem associated with (y*^'*)^\ (for given X^). Then, applying the standard argument of 
the Skorokhod problem, we get 



Lr'= sup f QHr.X'')dLl'^^ -Y,A , 



^ = l,...,iVi. 
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Hence, similar to Proposition [31 the same arguments as in the proof of Theorem 2 of [22] leads 
us to 

^ |Lf - Lf 'f < k sup 1^/'" - n*''f , < s < t (3.46) 

. , s<r<t ■ -. 

1=1 — — 1=1 

for some K > Q which depends only on V . 



Proposition 15. supE*^[ sup |XfT] < oo for all T > 0. 

M 0<t<T 

Proof. Take any /2 > and set tr = mi{t > ; |X*^| > R} and mf = E *'^[ sup |X 

0<r<min{t,rj;} 

From f l3.46p . we see that 



M|4l 



^ sup 



s<r<t 



1 ^1 



1 2 



'0 i=l 
N 



(3.47) 



hence [A6], Proposition [TU the Holder inequality, and the Burkholder-Davis-Gundy inequality 
imply 



R 



< cll + T^V l\''[l{r,>r} sup |6*^'X5,X^)|Vr 

JV /.t 1 

/ E"^[l{.«>.} sup |a^-^^(.,X^^)|Vr 

jit 0<s<r J 



< <! 1 + / m^dr \, t<T 



for some C > and > 0. Then we apply the Gronwall inequality to get = 0. Our 
assertion is now obtained by letting i? — )■ oo. ■ 

The inequality (13.461) and Proposition [15] immediately give the following proposition. 

Proposition 16. sup E^'^[(Iv^01 < ^ for all T > 0. 

M 

Proposition 17. The family of processes (X*^ , Z^' , L^'')m is tight on C^+^^+^i . 

Proof. By (I3.46p . (I3.47p . and a calculation similar to that in the proof of the above proposition, 
we have 



N Ni 



E^^niXr -X 



i=l i=l 

t 



M,i\Al 



<CVe'"[sup - y/^'f ] < C^(t - s) / (1 + E''[sup \X^\'])dr 

~^ s<r<t Js 0<v<r 



< C^(l + supE'^'[ sup \X^'\^]){t-s)\ 0<s<t<T 

M 0<r<T 



(3.48) 
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for some C, C'rp > 0. Similarly, 

E''[\Zt' - Zf\'] < C'^il + supE*'[ sup |Xf - sy (3.49) 

M 0<r<T 

for some > 0. The inequalities fl3.48p -( 13.49p . Proposition [TB| and Theorem 2.3 in [23] then 
give the assertion. ■ 

Proof of TheoremUl By Proposition [T71 we see that for any non-decreasing sequence {Mk)k 
there is a subsequence {Mki)i C {Mk)k and continuous processes {X,(p,Y, Z, L) on a certain 
probability space {Q, J-", P) such that 

(X^^i , ^^''"i , r*'^^; , , L^'"' ) — ^ (X, ip,Y,Z,L), / oo weakly. (3.50) 

We define {Ni)t by 

iV; = Z,^- / P\r,X)dr. (3.51) 

As in Proposition [H we get the weak convergence {X^^'^i , Z^'^i , L^^^; , iV^^*; ) — y (X, Z, L, X) 
by Proposition [T71 Thus, using fl3.38p - fl3.39p and the martingale representation theorem, we 
can find an X-dimensional (J-i)t-Brownian motion {Bt)t on a certain filtered space [Cl, J^, {J^t)t, P), 
which contains the original probability space {Q,J^,P), such that 

^t=y2 ^'Kr.X)dBi. (3.52) 

for i = 1, . . . , N. As in Proposition [T3|, we get 

F^o = xo + y —^dZ^, Y; = + Zl - f ^p^dZ\ ^ = 1, . . . , X. (3.53) 

Moreover, by fl3.43p . Proposition Uni and Theorem 2.2 of [T7], we get 

POO 

/ v^ic/Lj = 0, 2 = l,...,Xi. (3.54) 
Jo 

By f l330|) -f l334l) and Proposition [H we see that (X, L) is a solution of our SDER (EH])- 
Since [A 7] implies that the distribution of (X, ip, L) is uniquely determined, we get the weak 
convergence (X^, y?^, L^) — )■ (X, Lf, L) as M — oo. The proof of Theorem [1] can be 
completed using the arguments in step (vi) of the proof of Theorem 3 in [13]. ■ 
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